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A  METHOD  FOR  SOLVING  THE  INVERSE  PROBLEM  OF  TRANSONIC  CASCADE 
FLOW 


Cong  Zengjing 

Computing  Center,  Academia  Sinica 
Zou  Zixiang 

Institute  of  Engineering  Thermophysics,  Academia  Sinica 
Abstract 


This  paper  presents  a  method  for  solving  the  inverse  problem  of  transonic  cascade 
flow  in  (6,9)  coordinate  system.  Given  the  upstream  and  downstream  conditions  and 
the  velodry  distribution  on  cascade,  the  flow  field  it  obtained  by  solving  the  full  potential 
equation  with  velocity  as  the  unknown  with  the  line  relaxation  method.  Subsequently, 
the  field  is  transformed  back  into  the  physical  plane  and  we  obtain  the  profile  coordinates. 
Thickness  and  deflection  are  not  restricted,  since  the  full  potential  equation  is  used.  The 
numerical  result  of  the  inverse  problem  for  the  Hobson  cascade  flow  by  the  above  method 
shows  the  applicability  and  fast  convergence  of  our  method. 


I.  Preface 

In  the  last  ten  years,  many  numerical  methods  have  appeared 
for  the  analysis  of  transonic  cascade  flow  such  as  those  by 
Denton  [1,8],  MacCormack  [2],  FLIC  [3]  etc.  It  is  necessary  to 
obtain  the  velocity  distribution  on  the  blade  surface  for  given 
cascade  configuration  and  inlet  and  outlet  conditions.  However, 
the  calculation  time  of  this  method  is  too  long  and  moreover, 
when  there  are  small  blade  surface  configuration  corrections  near 
the  adjacent  velocity  line,  this  can  cause  tremendous  changes  of 
the  blade  surface  velocity.  Therefore,  to  obtain  a  cascade  with 
rational  velocity  distribution,  it  is  necessary  to  carry  out 
a  great  deal  of  molding  and  calculations. 


One  suitable  method  for  designing  a  cascade  with  low  loss 
[4]  begins  with  the  rationality  of  the  velocity  distribution 
and  as  regards  the  configuration  of  the  cascade,  for  the  velocity 
distribution  of  the  cascade,  the  experiences  have  been  voluminous 
(mainly  research  on  the  boundary  layer)  and  for  the  velocity 
distribution  and  the  inlet  and  outlet  conditions  of  a  given 
cascade,  they  find  the  configuration  of  the  blade  surface. 

This  then  composes  the  so-called  inverse  problem  of  cascade  flow. 

It  is  well  known  that  when  solving  the  mixed  potential  flow 
equation  using  the  potential  function  as  the  unknown  quantity, 
there  are  two  problems  worthy  of  attention.  One  is  finding  the 
hyperbolic  point  and  the  elliptic  point  positions  in  the  solution 
range  which  are  unknown  and  irregular.  If  the  calculation  mode 
used  is  not  suitable  this  will  produce  calculation  instability. 

The  other  problem  worthy  of  attention  is  the  use  of  the  physical 
boundary  condition  <p^= 0  is  not  convenient  especially  the 
problem  of  the  substance  surface's  form  being  too  complex. 

The  reason  is  the  grid  lines  included  in  the  substance  surface 
are  often  nonorthogonal . 

In  1971,  Murman  and  Cole  [5]  first  pointed  out  that  the 
"type  correlation  form"  caused  research  on  the  transonic 
potential  flow  problem  to  make  considerable  progress  yet  it  is 
only  limited  to  small  perturbation  equations.  When  the  flow 
direction  and  grid  line  direction  are  different,  the  negative 
artificial  viscosity  possibly  produced  in  the  supersonic  region 
causes  the  calculations  to  be  unstable.  In  1974,  Jameson  [6] 
pointed  out  that  the  "rotation  finite  difference  form"  eliminated 
the  limitation  of  the  small  perturbation.  However,  the  problem 
of  how  to  correctly  use  the  substance  surface's  boundary  condition 
#n*0  is  to  date  still  difficult  to  process.  When  the  wing  or 
cascade  thickness  is  very  thin,  4>  n«0  and  we  can  approximately 
substitute  with  ^^*0.  In  reality,  this  further  adds  to  the 
limitation  of  the  small  perturbation. 


This  paper  uses  velocity  q  as  the  unknown  quantity,  finds 
solution  in  the  (  V  ,  SP  )  coordinate  system  (similar  to  [4])  and 
readily  solves  the  above  two  problems.  Because  the  upper  and 
lower  sides  of  the  grid  here  are  flow  lines  themselves,  it  does 
not  require  structural  rotation  form.  Moreover,  because  the 
grid  is  rectangular  on  the  (  Y ,  )  calculation  plane,  this  is 

advantageous  for  structuring  the  implicit  form  advanced  along 
the  flow  line.  At  the  same  time,  the  orthogonality  of  the 
explicit  grid  naturally  satisfies  the  impermeable  conditions  of 
the  solid  wall  boundary. 

Symbols 

aT . :  critical  velocitv 
J-0 

s  :  arc  length  along  flow  line 
n  :  arc  length  along  potential  line 
'V'  :  flow  function 

:  inlet  and  outlet  angle 
r  :  specific  heat  ratio 
r  :  circulation 
w  :  velocity 
q  :  w/aLj 

P  :  density 
CP  :  potential  function 
&  :  local  flow  angle 
t  :  cascade  distance 

Lower  Symbols 

"1,2":  inlet  and  outlet  parameters 
"0"  :  stagnation  parameter 

II.  Fundamental  Equations 

Assuming  the  flow  satisfies  the  continuity  conditions,  flow 


function  ^  can  be  defined  as 

d<it  *■  pwdit 

Assuming  the  flow  satisfies  the  non-rotational  conditions, 
potential  function  sP  can  be  defined  as 

dip  “  u ids 

Because  the  dn  along  the  flow  line  is  zero,  when  the 
definition  of  is  known,  V”  is  constant.  The  ds  along  the 

velocity  potential  line  (vertical  to  the  flow  line)  is  zero  and 
when  the  definition  of  is  known,  is  constant. 

Naturally,  the  following  relational  formulas  exist: 

dx  co*  9ds  —  (co  *91  w)dq>, 

dy  —  niM  —  (ria0/ u>)d<p,  (2.1) 

dx  “  —  zln  Qdn  “•  ( —  nnd  I  pn/)d<& , 
dy  —  ,«mM*  *"  (eo *Bjp»')d<t>.  (22) 

Considering  the  two-dimensional,  non-rotational  and  constant 
compressible  flows 


G***.)  “  o. 

Os 

-5-  (•*<//)  -  o. 

On 


Further,  based  on  considerations  of  geometric  attributes 
(see  Fig.  1) ,  we  have 

_L  (2.5) 

*.  Ox  On 

1  d(Bs)  __  __  00, 

»,  On  Ox  * 

Based  on  formulas  (2.3),  (2.4),  (2.5)  and  (2.6),  we  obtain 
the  following  set  of  differential  equations: 
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-r-  (p«0*  +  6*  —  0, 
P* 


(2.7) 


f-  -  6,  -  0.  (2.8) 


From  (2.7)  and  (2.8)  ,  we  obtain  the  logarithmic 

differential  equation 

(lap)*,  +  (1 !«>)„  —  ((lap),)1  —  (lap)#(la»),  +  p*(lnp)^(la*-)# 

+  p’(W)**-0.  (2.9) 

We  introduce  the  isentropic  relationship 


(• _ i  \i/(r-i> 

‘"T+i*1)  • 


(2-10) 


Dimensionless  factors  p  ,  a  .  and  t  are  used  to  make 

O  Jj] 

formulas  (2.9)  and  (2.10)  dimensionless.  The  entire  potential 
equation  on  the  calculation  plane  (  y -  *F  coordinate  system) 
can  be  written  as 


*i(  !«»*)„  +  *i(l +  *i((ln*),]*  +  r«[(lnf)«]*  —  0. 

Formula  (2.10)  changes  into 

,-hri«‘)  •  a,2) 

In  the  formula,  c^=  (l-q2)/k ,  c2=l. 


(2.11) 


We  can  see  that  the  local  subsonic  and  supersonic  regions  of 
the  flow  separately  corresond  to  the  elliptical  and  hyperbolic 
regions. 
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Fig.  1  Geometric  relationship  under  (  V -  )  grid. 


III.  Boundary  Conditions 

The  boundary  conditions  of  the  solution  region  shown  in 
Fig.  2  are  given  as  follows: 

Inlet:  q^,  outlet:  q2,  (  $  2) ;  blade  surface:  qp,  qg. 

In  this,  ( /<5>  2 ^  is  the  re^erence  value. 

Letting  Q= (  p  2q2sin  /ff 2) / (  p  ^q^sin  /3 t  we  consider  that 

the  boundary  layer  causes  the  main  flow  region  of  the  wake  flow 
section  and  can  select  Q>1  so  that  the  wake  flow  section  is 
not  sealed. 

From  the  continuous  equations  and  isentropic  conditions, 
we  can  calculate 


ft 


train  {  \  Qqi  an  fi 


Fig-  2  Solution  region  and  boundary  conditions. 


Naturally,  when  Q  gives  fiz  ,  we  can  calculate  from  q  ,  $ 
and  q2- 

A  s^  and  A s2  on  the  pressure  boundary  flow  line  are 
calculated  as  follows: 


ix 


d, 

Min  A 


dl  — l  sin 


Circulation  V  along  the  cascade  profile  is  calculated 
based  on  strict  periodic  conditions 

T  —  —  ft  COS  A). 

If  under  design  conditions,  it  is  permissable  to  have  a 
supersonic  region.  The  necessity  to  thus  select  the  velocity 
distribution  causes  it  to  further  compress  the  isentropy. 

The  final  shock  wave  strength  will  be  weakened  and  thus  avoid 
the  separation  of  the  boundary  layer. 

It  is  only  necessary  that  the  boundary  layer  not  have 
separation  to  be  able  to  use  the  cascade  velocity  distribution 
with  the  largest  pressure  gradient. 


The  velocity  distribution  in  this  paper  is  given  along 
the  flow  line  isopotential. 

IV.  The  Discete  Technique 

When  there  is  discretization,  there  is  central  difference 
with  second-order  precision  along  the  direction.  However, 
along  the  SP  direction,  we  then  use  the  "type  correlation" 
form. 

The  elliptic  point: 

Cl®  “■  l  ~  21 &qa  + 

The  hyperbolic  point:  the  finite  difference  forms  of  first- 
order  and  second-order  precision  are  separately: 

(1®  ™  ~~  [lagfi  —  2 la $,■_!„  +■  ln^_2.j] 

(ln^)w|,j  «■  •"  * r  ■  [ 2 la 9,i  —  Slaf.-i.i  +  4lnf,_w  — 

In  order  to  avoid  processing  the  set  of  nonlinear  equations 
the  value  of  the  (n+1)  layer  only  appears  in  the  i  line,  on  the 
unknown  quantity  (Inq)  obtained  after  using  the  second-order 
difference  quotient. 

Finally,  we  obtain  a  set  of  finite  difference  equations 
which  assume  the  form  of  three  opposite  angles 

+  BjrOaf)r“  +  *«(  1®*).'-*”  -  ttV  (4.1) 

Coefficient  a!1^  etc.  are  calculated  as  follows 


Air  -  *\a'K&*y  -  *s? 

BIT  -  -2ri*7(A*)*  +  Hi-2c\‘' KA<py]  +  (1  +  H)[H. 

+  2(1  -  W,)M*7(A<p)’, 

£),*'  -  —  f  j*’ [ ( In  q.+uj  -  In  /(2  A?)  ]J 

—  fi*'[(ln?1„+,  —  In  /  C2  A<^»)  ]* 

—  H(ln9(+l„  ■+•  ln?,.ll>)f'i,7(A(p)a 

+  (1  -  H){W1[(21n^.w  -  lnf,_2./>i‘)/(A<p)J] 

•+•  (1  —  //i)[(51n^,_,.,  —  4ln^,_j.i  +  lflf, _».,)]}» 

w,  I1’  (1)S  q  <  *’ 

10,  q>  1. 

Key:  (1)  When. 

When  H^=l ,  formula  (1)  has  first-order  precision  (the  q 
direction)  and  when  H.^0,  formula  (4)  has  second-order  precision 

We  use  the  method  of  pursuit  to  solve  the  set  of  finite 
difference  equations  based  on  the  following  steps: 

For  a  fixed  i,  we  select  the  number  from  2  to  1-1  and 
carry  out 


1,  E!‘>  - -A? KB'.?  +  *'.?*£&’)» 

-  (D.r  -  +  c!?E}H), 

2.  ~  E  "  +  f;*'>  /  -  4  -  1- 

SB  (1) 

(2)  l 

IBA:  lvr..f,.P2.«A*fAl 
•  [»«*»*.)  ( 3 ) 


;  —  2-  •  —  l 


•2. 


r 


mtffSr.*****#®  (4) 
f5)  « 

In+Wltt  (lnf)lfZj 


CfcafifcJS  u«*“*  Oo#)**1)  (6) 

i 


L_ 

ll«.*B(2.1>!<2.2>*»,»]  (8) 

l 


Key:  (1)  Flow  chart;  (2)  Input;  (3)  Revised  velocity 
distribution;  (4)  Investigated  circulation  F  , 
if  not  suitable  then  change  to  2;  ^5)  Send; 

(6)  From  method  of  pursuit  use  Inq  to  calulate 

(Inq)n+^;  (7)  Use  q  to  calculate  0  along  middle 
flow  line  based  on  formula  (2.7)  and  calculate 
&  along  direction  based  on  formula  (2.8) ; 

(8)  Use  q,  &  to  calculate  x,  y  based  on  formulas 
(2.1)  and  (2.2);  (9)  Print  x,y  and  if  not  suitable 
then  change  to  2. 


After  the  calculation  of  the  entire  field,  we  then  further 
relax 


(la  —  <||(  la  ?))•*  +  (1  —  wXlaf)]^ , 

In  this,  relaxation  factor  Gj  6  [0,2]. 


The  total  calculation  procedure  is  shown  in  the  flow  chart. 


V.  Conclusion 


Firstly,  based  on  the  positive  problem  solved  in  Reference 
[7],  we  obtained  the  velocity  distribution  on  the  Hobson  cascade 
blade  surface  and  afterwards  used  the  method  proposed  in  this 
paper  to  sove  the  inverse  problem.  After  23  iterations,  the 
maximum  value  of  the  entire  field's  absolute  value  of  error 
8(\aq)  decreased  to  0.0002  and  we  obtained  the  profile  shown 
in  Fig.  3.  The  symmetry  accords  quite  well. 


Fig.  3  Calculate  Hobson  profile. 
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